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The spectrum of particles in supersymmetric Yang-Mills theory is expected to 
contain a spin 1/2 bound state of gluons and gluinos, the gluino-glue particle. We 
study the mass of this particle in softly broken supersymmetric Yang-Mills theory 
on a lattice by means of numerical simulations. The main focus is the estimation 
of finite size effects. We extrapolate the mass first to the infinite volume and 
then to the limit of a vanishing gluino mass. The results indicate that finite size 
effects are tolerable on lattices of moderate size, and that remaining deviations 
from supersymmetry are probably due to finite lattice spacing effects. 

1. Introduction 

Supersymmetric Yang-Mills theory (SYM) describes interacting gluons and their supersym- 
metric partners, the gluinos. As it is a gauge theory containing fermionic degrees of freedom, 
it is in this respect similar to QCD. An essential difference, however, is that gluinos are Ma- 
jorana fermions in the adjoint representation of the gauge group. In Minkowski space the 
(on-shell) Lagrangian of SYM is composed out of the gluon fields and the gluino field A, 
and reads 

l _ rn _ 

(1) 



tr 



i _ fji _ 

--Fu U F^ v + -X^D U X ^AA 

4 2 M 2 



where is the usual non-Abelian field strength and denotes the gauge covariant 
derivative in the adjoint representation. The fields A and are transformed into each 
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other by the supersymmetry transformation. The gluino mass term breaks supersymmetry 
softly. 

Some properties of SYM are expected to be similar to QCD [T]. It is asymptotically 
free and is assumed to show confinement. The determination of its low-energy properties, 
including the spectrum of particles, is a non-perturbative problem. The "physical" particles 
are bound states of gluons and gluinos, and if supersymmetry is unbroken, they would form 
supermultiplets . 

There are several motivations for the numerical investigation of SYM on the lattice. One 
of them is to understand the non-perturbative interaction of the gluinos in supersymmetric 
extensions of the standard model. There exist theoretical predictions for the low energy 
effective theory [2J E] that can be compared to the lattice results. Another motivation for 
the numerical investigations of SYM on the lattice is related to possible connections to 
ordinary QCD, as provided by the orientifold planar equivalence [3]. Results of previous 
investigations of SYM on the lattice by our collaboration can be found in [DJ . 

In recent years, SYM on the lattice has also been investigated with Ginsparg- Wilson 
fermions, in the domain wall [3 El E] as well as in the overlap formulation [ID]. For large 
lattice volumes and small lattice spacings these formulations require, however, a significantly 
larger amount of computing resources than the Wilson formulation. The gain of no need for 
tuning the position of the zero gluino mass point does not compensate by far the advantage 
of Wilson fermions. 

Our current studies of SYM are focussed on the bound states of gluons and gluinos, in 
particular on an exotic particle state in the spectrum that arises due to the fermions being 
in the adjoint representation. This particle is a spin 1/2 Majorana fermion and is called 
gluino-glue. It can be created by operators combining the field strength and the gluino field, 
the simplest example being 

6 si j = £<Vtr[i^A], (2) 

flU 

with a^u = \ [7^,7^]. Such a bound state containing a single fermion does not occur in 
QCD, but analogous particles exist in models similar to QCD with an arbitrary number of 
quark flavours in the adjoint representation. 

Since supersymmetry is generically broken in any non-trivial theory on the lattice [IT] 
it has to be ensured that it is restored in the continuum limit. A necessary condition 
for unbroken supersymmetry is the degeneracy of fermionic and bosonic masses. In a su- 
persymmetric theory the fermionic gluino-glue state must therefore be part of a multiplet 
containing also bosonic particles with the same mass. The behaviour of the gluino-glue 
particle provides an important signal for the supersymmetric limit of the theory. 

It is expected that in SYM a fine-tuning of the bare gluino mass parameter in the con- 
tinuum limit is enough to approach the symmetries of the continuum theory |12j . These 
symmetries include chiral symmetry and supersymmetry. The theoretical prediction of the 
existence of a supersymmetric chiral continuum limit needs to be confronted with the nu- 
merical lattice simulations. The chiral limit, the continuum limit, and the infinite volume 
limit can only be extrapolated from the results of these simulations. In QCD there are 
good estimates about the scales and parameters that are necessary to get reliable results, 
together with systematic estimates of the induced errors. In SYM the influence of para- 
meters like the finite volume and non-zero lattice spacing may be different. It is therefore 
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important to have a detailed understanding of their effects in order to be able to approach 
the supersymmetric limit. 

The previous results of our simulations have not yet shown the expected degeneracy of 
the fermionic and bosonic masses [5J E]. The obtained mass of the gluino-glue has been 
larger than the other masses of its lightest possible superpartners. However, the masses 
were obtained at a fixed lattice spacing and without a detailed analysis of the finite size 
effects. Our most recent investigations indicate that the influence of the finite lattice spacing 
is larger than expected. This provides a possible source of the supersymmetry breaking in 
the simulation. 

In this work we want to scrutinize the influence of the finite volume. If it is large, it 
provides another possible explanation for the gap between fermionic and bosonic masses. 
If it is smaller than expected then the lattice artifacts have to be reduced by performing 
simulations at smaller lattice spacings. 

In our studies the mass of the gluino-glue particle has typically been the one determined 
with the best accuracy. Therefore it is well suited to estimate the effects of the finite lattice 
size. In order to estimate the influence on supersymmetry breaking, the mass of the adjoint 
version of the rj meson (a-r/, a pseudoscalar bound state of gluinos) is also considered. 



2. Supersymmetric Yang-Mills theory on the lattice 

In our investigations we have employed the lattice formulation of SYM proposed by Curci 
and Veneziano [12]. The gauge field is represented by link variables U^(x) in the gauge 
group SU(-/V C ). The corresponding gauge action is the Wilson action built from the plaquette 
variables U p . The gluinos are described by Wilson fermions in the adjoint representation. 
In its basic form the complete lattice action reads 

S L = (l - ^fRetrcO + ± £ \ x (D w ) xy \ y , (3) 

p ^ c ' xy 

where D w is the Wilson-Dirac operator 

(J^w)x,a,a;y,b,/3 — $xy$a,b$a,fi 



ab3x—[i,y ■ (4) 

The hopping parameter k is related to the bare gluino mass via k = \/{2m g + 8). 

While the link variables U^{x) are in the fundamental representation of the gauge group, 
the variables V^x) in the Wilson-Dirac operator are the corresponding elements in the 
adjoint representation. In our present investigations the gauge group is SU(2), the funda- 
mental representation is a doublet and the adjoint one a triplet. In this case the adjoint 
gauge field is given 

by [Vf,{x)] ab = 2tr[C/t( x )r a C/ At (x)T 6 ], where T a are the generators of 
the gauge group, normalised such that 2 tr [T a T b ] = S ab . 

In our simulations the basic form of the lattice action has been modified with the following 
improvements: in order to reduce the lattice artifacts a tree-level Symanzik improved gauge 



3 



action has been used instead of the simple Wilson gauge action. One level of stout smearing 
has been applied to the link fields in the Wilson-Dirac operator to reduce lattice artifacts 
also in the fermionic part of the action. 

The lattice action explicitly breaks supersymmetry and in addition the chiral XJ(1)r sym- 
metry. However, to recover the continuum symmetries the necessary fine-tuning of super- 
symmetry and the U(1)r symmetry can be achieved through the same parameter, the bare 
gluino mass (i.e. k) |12j . The supersymmetric continuum limit coincides with the chiral 
continuum limit of the theory. 

An important similarity of SYM to Yang-Mills theory is confinement. The static potential 
between quark sources in the fundamental representation of the gauge group shows no signal 
of a string breaking. The low energy effective action is built out of bound states of the 
elementary fields. The glueballs of Yang- Mills theory are completed with the mesonic states 
(gluino-balls) and the fermionic gluino-glue particle. The proposed supersymmetric low 
energy effective actions are constructed from two multiplets of particles [2J [3] . The lighter 
multiplet consists of glueballs and a gluino-glue state. The heavier multiplet is built from 
the mesons and a gluino-glue state. As required by supersymmetry, each multiplet consists 
of a scalar, a pseudoscalar, and a fermionic spin 1/2 particle. 

On the lattice the ++ and the ^ glueball operators are constructed from the link 
variables. The mesonic particles are the a-/o represented by the operator A A and the a-r/ 
by A75A. Since the model contains only one species ("flavour") of gluinos, all possible mesons 
are "flavour diagonal". Consequently, their correlation functions contain disconnected parts 
in addition to the connected ones, as it is the case in QCD for flavour diagonal mesons. 

In addition to the above mesons, we consider the adjoint pion (a-7r), which is the pion in 
the corresponding theory with two Majorana fermions in the adjoint representation. The 
correlator of this particle is the connected contribution of the a-r/' correlator. The a-7r is 
not a physical particle in SYM, which only contains one Majorana fermion. However, it 
can be defined in a partially quenched setup, in which the model is supplemented by a 
second species of gluinos and the corresponding bosonic ghost gluinos, in the same way as 
for one-flavour QCD [13 . 

The reason for considering the a-n is the following. In the limit of a vanishing gluino 
mass, the bare parameter k has to be tuned to the critical value k c (chiral limit) that 
corresponds to the chiral theory in the continuum limit. The value of k c is most easily 
obtained from the dependence of the a-7r-mass on k. On the basis of arguments involving 
the OZI-approximation of SYM [2] , the adjoint pion mass is expected to vanish for a mass- 
less gluino. The a-7r yields a more precise signal for the tuning than the supersymmetric 
Ward identities. However, in previous studies we have checked that both signals lead to a 
consistent value of k c [S]. 

For both the a-/o and a-?/ mesons the correlators contain a significant contribution 
from disconnected diagrams (see Sec. |3]), especially on larger lattices and small adjoint pion 
masses. As in comparable cases in lattice QCD, this contribution leads to a bad signal to 
noise ratio in the correlators. Using the techniques detailed in Sec. El the gluino-glue mass 
can be obtained with an accuracy significantly better compared to the other particles of the 
spectrum (apart from the unphysical adjoint pion). 

In previous investigations the predicted multiplet structure of the particles has not been 
found in the mass spectrum [5] . The mass of the gluino-glue particle appeared to be heavier 
than the masses of the scalar and pseudoscalar states. A breaking of supersymmetry in the 
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theory on the lattice is induced by the finite lattice spacing and the finite lattice extent. 
The former is an unavoidable consequence of the discretisation; the latter is due to the 
anti-periodic (thermal) boundary conditions implemented for the gluinos. It is important 
to have an estimate of these effects. If they had an opposite sign, they could compensate 
each other and suggest a wrong supersymmetric point. Alternatively, they could sum up to 
the total supersymmetry breaking and it would be necessary to investigate both of them in 
order to disentangle their effects. 

In our numerical simulations we have applied a polynomial hybrid Monte Carlo (PHMC) 
algorithm. In a two step procedure, in the Metropolis step a better polynomial approx- 
imation is used than in the integrator of the molecular trajectory. The remaining error is 
compensated by a reweighting. Further details of the simulations algorithm can be found 
in jniig. 

The simulation at the smallest lattice spacings requires a large number of configurations 
to get a reasonable statistical error of the correlation functions. A further problem arises 
from the dependence of the adjoint pion mass on the lattice volume. This dependence is 
shown in Fig. [I] for a fixed value of k = 0.1490. The deviation of the smallest lattice from 
the infinite volume limit is most likely due to a larger influence of the excited states in the 
determination of the mass. 
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Figure 1: The mass of the adjoint pion (a-7r) at different lattice volumes for a fixed value 
of k = 0.1490. The fit (fit) is done using the ansatz Eq. [IT] and the reduced fit 
(red. fit) discards the smallest lattice volume, ro is the Sommer scale parameter. 



The adjoint pion mass gets smaller on the smaller lattices. Hence, for a fixed value of 
k the simulations are shifted towards the chiral limit when the lattice volume is reduced. 
In the vicinity of the chiral limit low eigenvalues of the Hermitian Wilson-Dirac operator 
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occur. The polynomial approximation in the PHMC algorithm has a larger error in this 
case and a compensation with correction factors is necessary. To obtain these factors we 
have calculated an approximation of the lowest eigenvalue on each of the configurations. 
Whenever this value has been below the threshold determined by the polynomial approxim- 
ation, the correct fermionic contribution of the 100 lowest eigenvalues has been calculated. 
From this contribution the correction factors of the reweighting step have been obtained. 
This reweighting leads to an additional increase of the statistical error on the smaller lattice 
volumes. A complete summary of the simulations is shown in Table 



3. Determination of masses 

The operator for the gluino-glue particle is represented by the lattice version of Eq. [5] 

°a = E<M P ^ A 1 ' ( 5 ) 

where the indices % and j stand for the spatial directions. In order to have the same 
properties with respect to parity transformation and time reversal, Fn V is represented by 
the anti-Hermitian part of the clover plaquette U^ c > [15] , 



*i = - (Uj$?) . (6) 

The clover plaquette is a combination of links in the fundamental representation, 

U$ = U^x)U v {x + ^)Ul{x + v)Ui{x) 

+Ul{x - u)Uf,(x - v)U v {x -u + n)Ul{x) 
+Uj l (x - n)Ul{x — fi — v)Ufj,(x - (j, — u)U u (x - v) 
+U u (x)Ul(x + u- fi)Ul(x - n)U^(x - fi) . 

The best signal is obtained from the contribution to the correlator that is proportional 
to the identity in Dirac space. The corresponding correlator is 

C gS {x -y ) = -| E E E E 

i,j,k,l x,y a,fi,p ab 

(a^tr[P l ,(x)T a ](D w - 1 ) x , a , /3 ,,, M tr[P fc/ (y)rV fc T)- (7) 
At large distances the correlator has the functional form 

C gS (t) » Csinh(?n(t - T/2)) , (8) 

where T is the temporal extent of the lattice. The mass m of the particle can be obtained 
by fitting the correlator to this function. The appropriate t-range for the fit, where t is 
large enough but still much smaller than T, can be found by plotting the effective mass 
m e ff(t). m e s(t) is the parameter m obtained from the correlator at t and t + 1 assuming 
Eq. [S] to be valid. In the t-range, where the influence of the excited states is small enough, 
m e ff(t) shows a plateau. In a more refined approach the correlated chi-square of the fit 



provides an indication that the assumption of the Eq. [5] is a good approximation in the 
considered region. We have applied a method that combines fit values of several fitting 
intervals and their correlated chi-square to get an estimation of the mass and the statistical 
and systematic error [TBI HI] • 

The gluino-glue correlator has been obtained using different smearing techniques. The 
link fields are smeared using APE smearing, the fermionic fields using Jacobi smearing. 
Without any smearing the links U in P^ should be the same fundamental links as in the 
gauge part of the action. The Wilson-Dirac operator D w , on the other hand, should contain 
the one-level stout smeared links in the adjoint representation, that are used in the fermionic 
part of the action. In addition to this approach, the one-level stout smearing has also been 
applied to the U fields. It has been checked that the smearing in the time-like direction 
included in this step has no influence on the correlation function at distances relevant for 
the mass determination. 




1 1 1 1 1 1 1 1 1 

012345678 

t 

Figure 2: The effective mass of the gluino-glue on a 16 3 x 36 lattice and k = 0.1492. The 
impact of different smearing parameters are compared. APE smearing is used on 
the the link fields in P, Jacobi smearing is applied on the fermionic source and 
sink. 

The effect of the smearing on the effective masses is shown in Fig. [5J The same Jacobi 
smearing is applied on the source and sink side of the correlator. The APE and the Jacobi 
smearing have a similar impact on the correlator. The excited states at small distances are 
considerably reduced by the smearing procedure. 
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The correlator for the a— rf boson contains a connected and a disconnected contribution, 



a 



Se-lf 



Ca-7r C a -^/ (Ji sc 

J3 Y.( tr [75(Dw)^75(Dw)^ 

X 

X^ tr [76(Dw)*i] tr [75(Dw)y,y 



(9) 



o 



Its connected part is the correlator of the adjoint pion (C a - n ). The disconnected part has 
been calculated using the stochastic estimator method jTS]- The statistical fluctuations in 
this part are usually larger than in the connected contribution. 
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Figure 3: The masses of the a-rj' and the adjoint pion (a— tt) at different lattice volumes for 
a fixed value of k = 0.1490. The fit (fit) is done using the ansatz Eq. [TTJ and the 
reduced fit (red. fit) discards the smallest lattice volume. 



It is instructive to compare the masses of the a-rj' and the a-ir, since their difference 
shows the influence of the disconnected part. In Fig. [3] the dependence of these masses on 
the lattice volume is shown. On the smallest lattice both masses deviate from the behaviour 
expected from the fits. This is presumably due to the larger contribution of excited states to 
the correlator, since only smaller values of t enter. A second observation is that at smaller 
lattice volumes both masses appear to approach each other. Thus at smaller volumes the 
influence of the disconnected part of the a-rj' correlator is relatively smaller compared to 
the connected one. 

At small t the contribution of the connected part is much larger than the disconnected 
part, whereas at intermediate distances and close to the chiral limit the disconnected part 
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yields a significant contribution. The larger statistical error of the a— 77' on the largest 
lattices is due to the larger statistical fluctuations of the disconnected part. The total error 
on a small lattice seems to be underestimated for the a— 7/ particle due to these systematic 
uncertainties. Smearing techniques or a larger temporal extent of the lattice could lead to 
better results. Because of these systematic uncertainties of the a-r/ mass and the larger 
statistical error, the best signal for the finite size effects is obtained from the gluino-glue 
mass. 



4. Finite size effects 

In quantum field theory in a finite volume the propagation and interactions of particles are 
different from those in an infinite volume. This can be understood as the modification of the 
polarized vacuum within a Compton wavelength around a particle due to the deformation 
by the finite boundaries of the box. The infinite volume mass mo of the particle is shifted 
due to the influence of the finite box size L, 

m(L) = m + Am(L) . (10) 

The mass shift Am(L) has been studied to all orders in perturbation theory in massive 
field theories in the continuum in A similar investigation in the framework of field 

theory on the lattice has been made in [2D]. The asymptotic behaviour was found to be 

Am(L) Rs CL" 1 exp (-am L) , (11) 

with parameters C and 1 < a < y/3/2. This behaviour is quite generic, since it does not 
depend on the specific form of the interactions. It also applies to the masses of stable bound 
states in a confining theory. For glueballs in lattice gauge theory the constants are 

3 V3 



c = ~w^ a = -' (12) 



where A is the three-glueball coupling constant. We fitted the dependence of the masses on 
the finite box size L by the general asymptotic relation Eq. [TTJ and obtained in this way 
an extrapolation to the infinite volume limit. To get more stable results, we have used a 
numerical fit to obtain the constants mo, C, and a. 

We have carried out simulations at several box sizes L, see Table [31 The temporal extent 
T of the lattice has been chosen to be about twice the spatial extent L. This combination 
has been chosen since it is the usual setup in all our simulations. 

The first estimation of finite size effects is done at a fixed value of the bare gluino mass 
defined by the hopping parameter k = 0.1490. To visualize the influence on the mass gap 
between bosonic and fermionic states, the gluino-glue mass is shown in Fig. 0] together with 
the mass of the bosonic a— rf meson. Our ansatz Eq. [TT] for the functional dependence of 
the mass on the lattice volume is valid only in the asymptotic region of large L and might 
fail for the smallest lattice sizes. To check its validity we have done a second fit (red. fit) 
that excludes the smallest lattice volume. The dimensionless scale 0.5L/ro corresponds to 
the length in femtometers if the Sommer parameter vq is set to ro = 0.5 fm as in QCD. Note 
that we have always taken the value of ro obtained by an extrapolation to the chiral limit 
of the data obtained on the 24 3 x 48 and 32 3 x 64 lattices. 
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Figure 4: The masses of the gluinoglue and the a— rj meson at different lattice volumes for 
a fixed value of k = 0.1490. The fit (fit) is done using the ansatz Eq. [Til an d the 
reduced fit (red. fit) discards the smallest lattice volume. 



The gluino-glue gets a positive mass shift at smaller volumes, whereas the a-r/' gets a 
negative one. This clearly shows that the finite lattice size enhances the supersymmetry 
breaking induced by lattice discretisation and nonvanishing gluino mass in our model: at 
smaller lattices an additional splitting of the bosonic and fermionic masses is introduced. 

In our previous simulations at larger box sizes the gluino-glue particle has been the 
heaviest of the measured low energy states. Since finite size effects are small in this case, 
this indicates that the mass shift induced by the discretisation artifacts is also positive. So, 
both finite size and finite lattice spacing effects add up to a total mass splitting between 
bosonic and fermionic states indicating the breaking of supersymmetry. 

In case of the gluino-glue, Eq. [TTJ seems to describe the volume dependence of the mass 
better than for the a-r/'. Also, the fit for the a-r/' fails for the smallest lattice volume. One 
possible reason for the different behaviour is the smearing of the gluino-glue that reduces 
the contribution of the excited states. This reduction becomes important for the smallest 
lattice sizes, where the masses are obtained from fits of the correlation function at smaller 
distances. Another reason are the systematic uncertainties in the mass determination of 
the a-r/' at smaller lattice sizes (see Sec. Nevertheless, neglecting the smallest lattice 
size, the fit seems to reproduce the general behaviour of the mass gap at different lattice 
volumes. 

Above a box size of about L = 1.2ro/0.5 (corresponding to 1.2 fm in QCD units) the 
statistical errors and the systematic errors of the finite size effects are of the same order and 
hence the finite size effects can be neglected at present. This is an important estimate for 
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the minimal lattice size required for the simulations. 

In the following we focus our discussion on the gluino-glue since the signal of the finite 
size effects is much clearer for this particle. Final results for masses in SYM should be 
obtained as extrapolations towards the infinite volume limit, the continuum limit, and the 
chiral limit. The continuum limit will be the subject of future investigations. The closer one 
gets to the chiral limit, the smaller are the masses of the particles and the ratio of the box 
size to the Compton wavelength. In other words, the finite size effects are expected to get 
larger in the vicinity of the chiral limit. Therefore, the proper procedure is an extrapolation 
to the infinite volume limit before the extrapolation of the chiral limit. 

In the first step we have extrapolated the values of the masses to the infinite volume limit 
but away from the chiral limit. Instead of a fixed value of the bare parameter k, we have 
fixed the squared mass of the adjoint pion, (rom a - T ) 2 , to five different values. The masses 
of the gluino-glue at these values have been obtained from a linear interpolation of the 
simulation results at each box size. These interpolations are shown in Fig. The masses 
at the three largest box sizes are almost indistinguishable in view of the current statistical 
errors. 

The results of the interpolation are shown in Fig. [6j The lines correspond to the fit of 
the dependence on the finite box size according to Eq. [TT1 The parameter tuq of the fit is 
the extrapolated infinite volume limit of the masses. The different parameters of the fits 
are summarised in Table [TJ 

The masses extrapolated to the infinite volume are shown in Fig. as a function of 
(?"o m a-7r) 2 - For comparison, the results of the largest lattice have been added to this plot. 

From the masses at different values of (rom^-^) 2 the chiral extrapolation can be obtained 
assuming a linear dependence. The linear regression of this dependence can be seen in Fig. [71 
The values obtained from linear extrapolations to the chiral limit are displayed in Table [2J 
The table contains the values at fixed lattice sizes as well as the results from the infinite 
volume extrapolations. The results at the largest two lattices appear to be consistent with 
those from the infinite volume limit. 



8 x 16 


12 x 24 


16 x 36 


20 x 40 


24 x 48 


32 x 64 


fit 


red. fit 


9.65(46) 


4.51(20) 


3.72(34) 


2.82(47) 


2.48(42) 


2.69(23) 


2.644(91) 


2.47(12) 



Table 2: The lower entries of this table are the results of chiral extrapolations of the gluino- 
glue mass at fixed lattice sizes. The last two entries are obtained by first performing 
the extrapolation to the infinite volume limit and then to the chiral limit. In the 
reduced fit (red. fit) the smallest lattice volume is ignored. All masses are in units 
of the Sommer scale. 
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Figure 5: The mass of the gluino-glue at different lattice volumes and different values of 
Ki clS cl function of the square of the adjoint pion mass in units of the Sommer 
scale. The lines are obtained from a linear regression of the points for each lattice 
volume. They are used as a linear interpolation. 
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Figure 6: A fit of the finite size effects at several fixed values of (rom a - x ) 2 . The fit is done 
using the ansatz Eq. fTTl The right panel shows the fits neglecting the smallest 
lattice volume. 



(roWa-Tr) 2 


TTlOr (fit) 


C (fit) 


a (fit) 


mor (red. fit) 


C (red. fit) 


a (red. fit) 


1 


2.85(13) 


33.6(71) 


0.68(11) 


2.74(19) 


10.8(84) 


0.39(26) 


2 


3.010(81) 


30.0(43) 


0.627(65) 


2.94(11) 


12.4(72) 


0.41(17) 


3 


3.165(48) 


25.0(28) 


0.555(46) 


3.108(60) 


14.0(47) 


0.426(92) 


4 


3.362(28) 


22.7(16) 


0.516(26) 


3.345(41) 


18.9(52) 


0.481(66) 


5 


3.543(56) 


22.3(33) 


0.510(54) 


3.578(70) 


35(23) 


0.60(15) 



Table 1: Fit parameters for the fits of the finite size effects shown in Fig. 
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Figure 7: The gluino-glue mass extrapolated to the infinite volume limit as a function of the 
squared mass of the adjoint pion in units of the Sommer scale. In addition the 
mass is shown at the largest lattice volume and different values of k. 



5. Conclusions 

In this work we have investigated the implications of the finite lattice extent in the simula- 
tions of SYM on a lattice. An exotic particle of SYM, the gluino-glue, turned out to provide 
the best signal for the finite volume dependence of the masses. A comparison with the bo- 
sonic a-r/ meson shows that the finite size effects on the masses have opposite signs for both 
particles: the gluino-glue gets a larger mass at smaller box size, whereas the mass of the 
a-r]' is reduced. This indicates that finite size effects enhance the supersymmetry breaking 
induced by lattice discretisation and nonvanishing gluino mass, and increase the unexpected 
gap between the bosonic and fermionic states seen in previous lattice results. The influence 
of the finite box size decreases, however, quite rapidly. The differences between the three 
largest lattices in our simulations are nearly negligible. The physical reason for this effect 
is the large mass of the lightest particle in this theory. The adjoint pion, which is light 
and becomes massless in the chiral limit, is not a physical particle and is only defined in 
a partially quenched framework. The lightest physical particles have significantly larger 
masses, even in the chiral limit. Therefore their Compton wavelengths are small compared 
to the sizes of our largest lattices. 

The aim of these investigations was to obtain estimates of the finite size effects and their 
implications on the setup for future simulations of SYM. Another item would be to consider 
the effect of changing the fermionic temporal boundary conditions from antiperiodic to 
periodic. This could help to disentangle the effects of the finite spatial volume and the 
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finite temporal extent. In our present work we have chosen a ration of L over T similar to 
our final simulations for the determination of the complete spectrum on larger lattices. 

The findings have implications for our further investigations of SYM. We have found that 
finite size effects are small for L > 1.2ro/0.5 (corresponding to 1.2 fm in QCD units), and 
that the simulations can be efficiently done at present lattice spacings on medium lattice 
sizes, for instance, 24 3 x 48. In order to obtain reliable results relevant to the continuum 
limit a sufficiently large statistics and a small lattice spacing appear to be crucial. 
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A. Summary of the simulations and results 
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Table 3: Summary of our simulations at /3 = 1.75. The value of r^/a obtained from the 
extrapolation of the 24 3 x 48 and 32 3 x 64 lattice data to the chiral limit is 
r o/ a = 8.89(10). iV con f is the number of configurations in the gluino-glue mass 
measurement. iV corr is the total number of configurations with a reweighting factor 
different from one. 
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